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The temperature evolution of the chiral condensates in a gas made of pious, kaous and etas is studied within 
the framework of SU(2) and SU(3) Chiral Perturbation Theory. We describe the temperature dependence of the 
quark condensates by using the meson meson scattering phase shifts in a second order virial expansion. We find 
that the SU(3) formalism yields an extrapolated melting temperature for the non-strange condensates which is 
lower by about 20-30 MeV than within SU(2). In addition our results show that the strange condensate melting 
is slower than that of the non-strange, due to the different strange and non-strange quark masses. 



1. Introduction 

We review here our recent study of the thermal 
evolution of the chiral strange and non-strange 
condensates [1]. Our motivation is the study of 
the the phase diagram of QCD at zero baryonic 
densities, which is a state that is expected to oc- 
cur in the central rapidity regions of Relativistic 
Heavy Ion Collisions. Our purpose is to obtain a 
model independent description of the condensate 
and its thermal evolution. To that aim we will 
turn to the virial expansion and to Chiral Pertur- 
bation Theory. 

The virial expansion is a simple and successful 
approach to describe many thermodynamic fea- 
tures of dilute gases made of interacting pious [2] 
and other hadrons [3]. Let us note that for most 
thermodynamic properties it is enough to know 
the low energy scattering phase shifts of the par- 
ticles, which, in principle, could be taken from ex- 
periment, avoiding any model dependence. How- 
ever, we are interested in the chiral symmetry 
restoration, whose order parameters are the quark 
condensates, which are defined as derivatives of 
the pressure with respect to the quark masses. 
These derivatives cannot be measured experimen- 
tally and therefore a model independent theoreti- 
cal description of the scattering amplitudes mass 
dependence is required. That is why we need the 
second ingredient of our approach: Chiral Per- 



turbation Theory (ChPT) [4-6], which provides a 
model-independent and systematic description of 
low energy hadronic interactions. ChPT follows 
from the identification of the pions, kaons and the 
eta as the Goldstone Bosons of the QCD spon- 
taneous chiral symmetry breaking (pseudo Gold- 
stone bosons indeed, due to the small light-quark 
masses). Then, the ChPT Lagrangian is built 
as the most general derivative and mass expan- 
sion, over 47ri^ ~ 1.2 GeV, (the symmetry break- 
ing scale) respecting the symmetry constraints. 
At one loop any calculation can be renormalized 
in terms of just an small set of parameters, Lk{fJ.), 
Hk (fJ.) , (/i being the renormalization scale) which 
can be determined from a few experiments and 
used for further predictions. 

In particular, within SU(2) ChPT the lead- 
ing coefficients of the low temperature expan- 
sion for the pressure and the non-strange quark 
condensate have been calculated for a hadron 
gas whose only interacting hadrons were the pi- 
ons [7]. Adding other more massive particles 
(hke kaons, etas, rhos, protons, etc..) as free 
states, it was possible to obtain an estimate, 
T ~ 190 - 200 MeV, of the phase transition criti- 
cal temperature. It was also shown that the per- 
turbative calculation of the pressure was analo- 
gous to the second order virial expansion when 
the interacting part of the second virial coefficient 



is obtained from the one loop tttt ChPT scattering 
lengths. Other works have studied the applicabil- 
ity of the virial expansion with a pionic chemical 
potential [8], or the critical temperature in gen- 
eralized scenarios [9]. 

The interest of extending the previously com- 
mented work to SU(3) is that the SU(2) approach 
is limited by the absence of other interacting 
particles, like kaons or etas, whose densities at 
T ~ 150 MeV are significant [7]. In addition, 
let us recall that the chiral phase transition can 
be different for the SU(2) and SU(3) cases, and 
that several QCD inequalities and lattice results 
suggest an stronger chiral condensate tempera- 
ture suppression with an increasing number of 
light flavors [10]. Intuitively this is due to the 
fact that the existence of other states (strange 
quarks, or strange mesons) makes it easier to cre- 
ate entropy, that is, disorder, and therefore it is 
easier to melt the ordered state, that is, the con- 
densates. This effect has been observed in lattice 
calculations [11], but only in the chiral limit, and 
it seems to lower the chiral critical temperature 
down by roughly 20 MeV. Note, however, that 
none of these results have been obtained from 
the hadronic phase and realistic masses, which 
is what we have achieved with our approach. 

In addition, within SU(2) it is not possible to 
study the (ss) condensate. This could be interest- 
ing because, in the chiral phase diagram, quark 
masses play the same role as magnetic fields in 
ferromagnets: Intuitively, we need a higher tem- 
perature to disorder a ferromagnet when there 
is a magnetic field aligned along the direction 
of the magnetization. Analogously, it has been 
found that the SU(2) chiral condensate melts at 
a lower temperature in the massless limit [7]. In 
particular, we have obtained the temperature de- 
pendence of (ss) in the hadronic phase. Since 
the strange quark mass rus is much larger than 
m„,md, we have indeed found a sizable "ferro- 
magnetic" effect which translates into a slower 
thermal (ss) evolution. This effect had already 
been studied in the large Nc limit with the use of 
QCD-motivatod effective Lagrangian [12]. 

Thus, the second order relativistic virial expan- 
sion of the pressure for an inhomogeneous gas 



made of three species: i = tt, /sT, ry is [13,2]: 

i i \ jjti J 

where P = 1/T and £,i = exp(— /3mi). Expanding 
up to the second order in means that we con- 
sider only binary interactions. For a free boson 
gas, = for i ^ j, whereas: 
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The degeneracy is gi = 3,4,1 for Tr,K,r], respec- 
tively. The interactions appear through [13,7]: 
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where Ki is the first modified Bessel function and: 
A^^-= ^(2/+l)(2J+l)<5)^;,_^(i?), (4) 

I,J,S 

Sfjg being the ij ij phase shifts (chosen so 
that (5 = at threshold) of the elastic scatter- 
ing of a state ij with quantum numbers /, J, S 
(J being the total angular momentum and S the 
strangeness). 

The virial expansion breaks around T ~ 200 — 
250 MeV [8]. Fortunately, the physics we are 
interested in occurs already at T < 250 MeV, 
and > — ^7], and it is enough to consider 
ij = TTTT, ttK and nr] in Bij. 

At this point we recall that the quark mass ap- 
peals in the QCD Lagrangian, and therefore in 
the partition function, as mq^qaqa, where qa = 
u,d,s. Thus, in order to obtain the condensate 
we simply have to differentiate the partition func- 
tion with respect to niq^ . At finite temperature, 
the partition function is substituted by the free 
energy density z, so that [7] 
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Note that we have separated the T = part from 
the temperature dependent part, which is nothing 



but the pressure P = eo — z, cq being the vacuum 
energy density. Of course, at T = we recover 
(qaqa) = {0\qaqa\0) = dea/dniq^. 

Once more, we emphasize that in contrast with 
most thermodynamic quantities, for the chiral 
condensate it is not enough to know the S{E), 
but we also need their dependence with the quark 
masses as well as a value for the T = vacuum ex- 
pectation value. Wc will obtain this information 
from ChPT. Nevertheless, ChPT does not deal 
with quarks, and thus we first have to rewrite the 
condensate, eq.(5), in terms of meson masses: 

( ifcqc. Qp\ 

where, as before, i = 'K,K, t], and we have defined: 
cr"=-F'^mqaqa\0)-\ (7) 

In the isospin limit, when to„ = mj,, the u and d 
condensates are equal, and we define (OlgglO) = 
{0\uu + dd\0). It is tedious but straightforward 
to obtain the c parameters above, whose explicit 
expressions can be found in [1]. The only relevant 
comment is that the c coefficients depend on the 
chiral parameters Lk, for k — 4. ..8, and i?2. 

In Table I, we show two Lj. determinations 
from meson data, but it makes little difference 
to use other parameter sets. For H2 wc will 
use H^iMp) = (-3.4 ± 1.1)10-3, obtained as ex- 
plained in [14] but using a more recent estima- 
tion of (0|ss|0)/(0|gg|0) = 0.75 ± 0.12 [15]. For 
instance, the c parameters obtained when using 
the Lk in the first column of Table I are: 

— U.y_o.4, — U.0_o.7' S — U.4_g 7, 

The above c|* is in good agreement with the 
SU(2) estimates: 0.85 and 0.90 ±0.05 [7]. 

Finally, wc have to introduce the meson-meson 
interactions in the second order virial coefficients 
through the phase shifts in eq.(4). These are 
nothing but the complex phases of the meson- 
meson amplitudes, once they are projected in par- 
tial waves of definite isospin and angular momen- 
tum. The one loop ChPT amplitudes have been 





Refs.[5,21] 


lAM [17] 


L\ 


0.4 ±0.3 


0.561 ±0.008 (±0.10) 


T r 
^2 


1.35 ±0.3 


1.21 ±0.001 (±0.10) 




-3.5 ± 1.1 


-2.79 ±0.02 (±0.12) 


LI 


-0.3 ±0.5 


-0.36 ±0.02 (±0.17) 


rr 
^5 


1.4 ±0.5 


1.4 ±0.02 (±0.5) 


LI 


-0.2 ±0.3 


0.07 ±0.03 (±0.08) 


Lr 


-0.4 ±0.2 


-0.44 ±0.003 (±0.15) 


LI 


0.9 ±0.3 


0.78 ±0.02 (±0.18) 



Table 1 



Different sets of chiral parameters xlO^. In the 
first column L^, L2, are taken from [21] and 
the rest from [5] (L4 and are estimated from 
the Zweig rule). The last column is the lAM fit 
[17] to meson-meson scattering up to 1.2 GeV. 



given in [16,17]. As already commented within 
SU(2) ChPT it has been shown [7] that using just 
the amplitudes at threshold (scattering lengths) 
in the virial expansion is equivalent to expanding 
the partition function to third order in T/F or 
T/m-r^. This approach yields a fairly good repre- 
sentation of the pion gas thermodynamics at low 
temperatures, T < 150 MeV [7]. Let us recall 
that ChPT provides a good low energy descrip- 
tion {E < 500 MeV) of the meson-meson ampli- 
tudes [16]. 

The virial coefficients and dP/drui have been 
calculated numerically. For simplicity, in the 
figures, we have represented the chiral con- 
densate over its vaciium expectation value, 
(9a<?a)/(0lga9a|0)) SO that all of them are nor- 
malized to 1 at T = 0. As a matter of fact, since 
there is always an small explicit chiral symmetry 
breaking due to the quark masses, the condensate 
should only vanish completely in the T — > 00 limit 
(following with our analogy between the quark 
masses and the magnetic field, our ferromagnet 
above the Curie point becomes paramagnetic, but 
as if it was in the presence of a magnetic field, 
which still produces some magnetization.). Of 
course, with a virial expansion truncated at sec- 
ond order we cannot generate such an analytic 
behavior, and our curves become negative above 
some T, where the approach becomes clearly un- 
reliable. Nevertheless, since the u and d masses 
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Figure 1. Condensate thermal evolution using one-loop 

ChPT amplitudes. Although they should not really vanish 
at a finite temperature, they are extrapolated down to zero 
only for reference. 



are so small compared with the size of the non- 
strange condensate at T=0 (less than 3%), it is a 
fairly good approximation to speak about a melt- 
ing temperature in that case. For the (ss) con- 
densate, we will only consider an extrapolated 
melting temperature to case the comparison with 
the non-strange case, but that number does not 
have an actual physical meaning. 

In Figure 1 we show the results of using the 
one-loop ChPT phase shifts obtained using the 
central values of the parameters in the first col- 
umn of Table I. For illustration, we have sepa- 
rated the different contributions. In particTilar, 
the thin continuous line represents the interact- 
ing pion gas in SU(2), (see [6] for the conversion of 
SU(3) into SU(2) parameters), whereas the thin- 
dashed line is the corresponding one-loop SU(3) 
result also for a pion gas. We have checked that 
the tiny difference between them comes only from 
the pure 0(p'*) phase shift contribution, ( the am- 
plitude, once the parameters arc translated, is the 
same for SU(2) and SU(3) only up to s/M| or 
s/M^ terms coming from the expansion of kaon 
or eta loops, which, of course, are not present in 
SU(2) [6]). This amounts to a 4 MeV decrease of 
the pion gas extrapolated melting temperature: 
T^"^^ = 231 MeV for SU(2) and Tif«> = 227 MeV 
for SU(3). Next, the thin-dotted line is obtained 
by adding free kaons and etas to the pion gas, 
which lowers the extrapolated melting tempera- 
ture, by roughly 5-6 MeV down to 221 MeV. 



Our SU(3) result for a. Tr,K,ri gas, where also 
the kaons and etas interact with each other and 
the pions, is the thick continuous line, where we 
see that the additional decrease, basically due to 
the TrK,-jTr] interactions, amounts to 10-11 MeV. 
On a first glance, it may seem surprising that the 
ttK, ttt] interaction effect comes out comparable 
or larger than that of free kaons and etas, since it 
is thermally suppressed by another exp(— m^/T) 
factor. Note, however, that the exp(— m,r/T) sup- 
pression only amounts to a factor of 6, 4, 2.5 at 
T= 80,100,150 MeV, respectively. In contrast, 
the ttK and tttj interactions depend strongly on 
m^, which is much more sensitive to m than 
rriK or m^, which carry the only to dependence 
of the free K and t] terms. In particular, it 
is easy to see that dm^/drh ~ 2dm%/dm and 
drriT^ / dm ~ Q drriK / dm, and this "temperature 
independent enhancement" competes with the 
thermal suppression, making the -kK and ttt/ in- 
teraction effect comparable to the free one, al- 
ready at T = 80 — 100 MeV and even larger if 
T > 140 MeV. 

Putting all the pieces together, the Tm"^^ de- 
crease in SU(3) is roughly 20 MeV, in agreement 
with the chiral limit lattice results quoted in [11]: 
Tc = 173 MeV for SU(2) and = 154 MeV, for 
SU(3). The effect of using the real masses would 
be to increase both temperatures. 

Furthermore, it can be noticed that (ss) melts 
much slower than {qq), since TOs ^ to. Indeed, 
there is still 70% left of the (ss) condensate at the 
{qq) melting point. 

Let us recall that both effects are already siz- 
able at low temperatures T ~ 100 MeV, where we 
can still trust pure ChPT. As commented above, 
only for illustration and to ease the comparison 
between curves and with previous works, we have 
extrapolated the condensates down to zero. 

At this point we can also consider the uncer- 
tainties in the chiral parameters in Table I. By 
performing a Montecarlo gaussian sampling of the 
parameters within their errors, we can obtain a 
rather conservative estimate of our uncertainties 
(since the sampling is done assuming uncorrelated 
errors). We thus find: 

Tjf«> = 21i;^^MeV, T^'^i = 291^^^ MeV 



Note that the (ss) extrapolation is beyond the 
reliability region of the second order virial expan- 
sion, is just illustrative, and really only vanishes 
in the T oo limit. Let us also remark that the 
extrapolated melting temperatures are strongly 
correlated, so that their difference has less un- 
certainty than simply adding in quadrature their 
respective errors. In particular, we find 

yj«g>,SC/(2) _ r^(^<l),SUi3) ^ 211^4 MeV 

We remark once more that the (ss) melting tem- 
perature is, of course, only a crude extrapolation 
for illustrative purposes and to ease the compar- 
ison with the non-strange condensate evolution. 

Furthermore, we have estimated the effect of 
other, more massive, hadrons. In the SU(2) case 
[7,8] these states also included the kaons and etas, 
although assuming a rather large uncertainty in 
dMi/dmq^. Ah in aU this effect reduced T^'^'' by 
approximately 10-20 MeV. In the SU(3) case we 
have indeed calculated explicitly the kaon and eta 
m dependence in the c'"?" parameters, thus re- 
ducing considerably the above uncertainty. Since 
the kaons and eta are treated explicitly, the other 
massive hadrons are heavier than m^, have a 
very low density, and their main contribution to 
the pressure comes from the first virial coeffi- 
cient, i.e. the free gas. The only uncertainty 
is on dMh/dniq^, conservatively estimated to lie 
within the number of valence quarks Nq^ and 
2Nq^ . Thus, we have found that their contribu- 
tion decreases T^f^^ by 7-12 MeV, respectively. 

Finally, and in order to estimate the effects 
of higher energies we can extend ChPT up to 
E ~ 1.2 GeV by means of unitarization models 
[18,19,17]. These techniques resum the ChPT 
series respecting unitarity but also the low en- 
ergy expansion, including the mass dependency. 
In particular, it has been shown that the cou- 
pled channel Inverse amplitude Method (lAM) 
provides a remarkable and accurate description 
of the complete meson-meson interactions below 
1.2 GeV, generating dynamically six resonances 
from the one-loop ChPT expansion and a set of 
fitted Lk compatible with other previous deter- 
minations. This approach can be extended sys- 



tematically to higher orders (for SU(2) case, see 
[20]). 

Hence, in Figure 2 we show the results of using 
the one-loop coupled channel lAM fitted phase 
shifts. Its corresponding Li parameters are given 
in the last column of Table I, with two errors, 
the first, very small, is purely statistical, and the 
second covers the uncertainty in the parameters 
depending on what systematic error is assumed 
for the experimental data. Let us remark that, 
although the Lk are highly correlated, the sec- 
ond, larger error, ignores completely these corre- 
lations, and should be considered as a very con- 
servative range. The continuous line corresponds 
to the central values, and the dark shaded areas 
cover the one standard deviation uncertainty due 
to the small errors in the parameters. These areas 
have been obtained from a Montecarlo Gaussian 
sampling. The conservative ranges are covered by 
the light gray areas. We now find: 

Ti«> = 2041:? (i|)MeV,rif'*> = 304^^9 {^l^)Me\ 

where the errors in parenthesis are the conser- 
vative errors which should be interpreted as un- 
certainty ranges better than as standard devia- 
tions. Note the excellent agreement with stan- 
dard ChPT. The magnitude of the different con- 
tributions is roughly the same, although the ttK 
and TTTj interactions this time lower TiP^ by 17 
MeV, and the free kaons and etas by roughly 
10 MeV. Again other massive states apart from 
kaons and etas would lower riP^ by 6-10 MeV, 
and tI^^^ by 22-32 MeV. In addition, we find 

Ar„ = rif«> - rif«> = lootl^ ( so ) MeV, 

The lAM results show that higher energy effects 
do not affect very much our conclusions and that 
the ChPT extrapolation from low temperatures 
seems robust. 

To summarize, we have studied the SU(2) and 
SU(3) temperature evolution of the chiral conden- 
sates in the hadronic phase. The thermodynamics 
of the meson gas has been obtained from the virial 
expansion and Chiral Perturbation Theory. Our 
results clearly show a significant decrease, about 
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Figure 2 Temperature evolution of cliiral condensates us- 
ing the unitarizcd ChPT amplitudes. The shaded areas 
cover the uncertainties in diflferent sets of chiral param- 
eters. Although they should not really vanish they are 
extrapolated down to zero only to ease their comparison 
and for reference. 

20-30 MeV, of the non-strange condensate melt- 
ing temperature, from the SU(2) to the SU(3) 
case, similar to lattice results. Of these, about 6 
MeV had already been explained with free kaons 
and etas, but the rest are mainly due to -nK and 
Trry interactions. We have also estimated the effect 
of heavier hadrons and of the third order virial 
coefBcient. All in all we find 

Tjf«> = 20lt?? MeV, 

In addition, our results show an slower temper- 
ature evolution of the strange condensate, shifted 
by about 80 MeV with respect to the non-strange, 
due to the different quark masses. More quanti- 
tatively, {qq) does not show a sizable melting up 
to T ~ 150 MeV and it still remains about 80% 
when the non-strange is basically zero. 

Both effects are clearly visible already at low 
temperatures. However, we have checked that 
their size is completely similar when using uni- 
tarized models. 

These techniques should be easily extended to 
Heavy Baryon Chiral Perturbation Theory, in 
order to study the condensates with non-zero 
baryon density. In general the whole approach 
could be used with any effective Lagrangian for- 
malism, in particular to study other QCD phase 
transitions like those of the color superconducting 
phases. 
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